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We demonstrate the existence of an exact invariant solution to the Navier-Stokes equa- 
tions for the asymptotic suction boundary layer. The identified periodic orbit with a 
very long period of several thousand advective time units is found as a local dynamical 
attractor embedded in the stability boundary between laminar and turbulent dynamics. 
Its dynamics captures both the interplay of downstream oriented vortex pairs and streaks 
observed in numerous shear flows as well as the energetic bursting that is characteristic 
for boundary layers. By embedding the flow into a family of flows that interpolates be- 
tween plane Couette flow and the boundary layer we demonstrate that the periodic orbit 
emerges in a saddle-node infinite-period (SNIPER) bifurcation of two symmetry-related 
travelling wave solutions of plane Couette flow. Physically, the long period is due to a 
slow streak instability which leads to a violent breakup of a streak associated with the 
bursting and the reformation of the streak at a different spanwise location. We show that 
the orbit is structurally stable when varying both the Reynolds number and the domain 
size. 
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1. Introduction 

Recent progress in understanding transitional turbulence and the transition itself is 
based on the application of concepts from dynamic systems theory. The image that 
emerges is centered around invariant solutions of the underlying Navier-Stokes equations 
such as fixed points, travelling waves and periodic orbits. Together with their entangled 
stable and unstable manifolds these solutions form the backbone of a chaotic saddle that 
supports the observed turbulent chaotic dynamics: a trajectory approaches the invariant 
solutions along their stable manifolds and leaves again along the unstable ones; it consists 
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of a chaotic sequence of close visits to the invariant solutions. Support for this concept 
comes in particular from numerical observations in confined geometries, including plane 



Couette flow (Nagata 1990 Kawahara fc Kida 2001 Viswanath 2007 Gibson et al. 2008 a 



Halcrow et aL|2009|, Taylor Couette flow (|Faisst k Eckhardt|2000[) and most prominently 
pipe flow (|Faisst fc Eckhai^[2003{ [Wedin fc Kerswell||2004| |Eckhardt et ctZ.||2008[ |Willis| 
et al. 2010), where flow structures very similar to the numerically calculated ones have 



been observed experimentally by Hof et al. ( 2004 1 . 



These flows are known to show a transition to turbulence for flow rates where the 
laminar profile is stable against infinitesimal perturbations ( Grossmann|2000" ) ; finite am- 
plitude perturbations are required to trigger the transition. In these systems, the study of 
dynamically attracting objects within the boundary that separates laminar from turbu- 
lent dynamics - the so called edge states - has proven very fruitful as these states seem to 
guide the transition to turbulence ( Toh fc Itano|20 03 ; Skufc a et aL|2006 Schneider et al. 



2007 2008). In some confined flow geometries, edge states haven been very well charac- 



terized; they can be invariant solutions like fixed points in plane Couette flow, travelling 



waves in symmetry reduced pipe flow or extended plane Couette flow (Schneider et al. 
2010) or be more complex and even chaotic in full pipe-flow ( Schneider et a/.||2007 ). 



The boundary layers that form around bodies in open flows show a somewhat richer 
dynamics than confined shear flows. In addition to near-wall downstream vortices and 
streaks, energetic bursting events occur where low-speed fluid is ejected from the wall into 
the free-stream. These bursts are responsible for the majority of turbulence production 
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to which extent the invariant solutions can capture both the near-wall vortex structures 
and the bursting. 

The transition to turbulence in boundary layers that form around bodies in open 
flows shares with the internal flows just described the fact that turbulence is observed 
in parameter ranges where the laminar flow is still stable. Techniques of edge state 
tracking, initially developed for internal flows, have been adapted and applied to the 



Blasius boundary layer by Cherubini et al. (2011) and Duguet et al. (2012). However, 



limitations on the domain size prevented a full characterization of the invariant solutions. 
More generally, the extension of the concepts developed for internal flows to external flows 
is interesting because of the rich phenomenology seen in external flows. For instance, in 
addition to the near-wall downstream vortices and streaks, which seem similar to the 
structures found in internal flows, energetic bursting events occur where low-speed fluid 
is ejected from the wall into the free-stream. These bursts are responsible for the majority 
of turbulence production in the boundary layer ( Kline et «Z.||1967 Robinson] 1991 1. The 
extend to which the invariant solutions can capture both the near-wall vortex structures 
and the bursting in boundary layers remains open. 

As a step towards open flows, we here investigate the turbulence transition in the 
asymptotic suction boundary layer (ASBL), a boundary layer with a mean normal flow 
(Fransson 2001 ) that shows many features of both internal and external flows. It is ideally 



suited for an investigation of invariant solutions and edge states in boundary layers, 
as it has a translationally invariant base flow. Furthermore, we gain the possibility to 
connect the system to plane Couette flow by tracking a homotopy parameter, which will 
allow us to fully characterize and explain the dynamics of the edge state. Finally, the 
flow connects to engineering applications in which suction is a well studied concept for 



turbulence control, see for example Joslin ( 1998 ) 



This paper is organized as follows: In |2]we describe the flow, its numerical implemen- 
tation and the algorithms used to find edge states. In Sj3]we discuss the flow fields for 
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Re — 400, followed by a discussion of the Reynolds number and aspect ratio variations 
in !j4j Concluding remarks are given in <|5] 



2. System, numerical scheme and algorithms 

2.1. The asymptotic suction boundary layer 

The asymptotic suction boundary layer (see figure [lj describes a situation where fluid 
flows with velocity Uoo over a flat plate into which it is sucked with a perpendicular 
velocity Vs- We work in a Cartesian coordinate system, with x the downstream direction; 
y and z are the wall-normal and spanwise directions, respectively. The equations of motion 
for an incompressible fluid are the Navier-Stokes equation 

(d t + u • V)u = -Vp + (/Au (2.1) 

together with the continuity equation 

V-u = 0. (2.2) 

The no-slip boundary conditions at the wall, the downstream flow and a constant, 
homogeneous suction Vs give for the stationary laminar solution to the full Navier-Stokes 
equation an exponential profile, 

u = (Uoo(l - c-v/t), -V s , o). (2.3) 

The length 5 = v/Vs equals the displacement thickness, calculated from J^°(l— u(y)/U 00 )dy, 
and is used as length-scale. Velocities are measured in units of and time in units of 
t u = 5/Uoq. The Reynolds number for the ASBL is defined as 

Re = ^ = (2.4) 
v V s 



The laminar profile (equation 2.5) is linearly stable for Reynolds numbers up to a 
critical value of Re c — 54382 ( Hocking|1975 ). This is a factor of about 100 more than for 



a boundary layer without suction Schlichting ( 2004 ) . As calculations in this work do not 



exceed Reynolds numbers beyond 750, the laminar flow is stable for all our considerations. 

For the numerical treatment, the unbounded wall-normal direction is cut off at a height 
h, where a second porous plate that moves with constant velocity is introduced. At 
both walls we enforce no-slip boundary conditions for the x- and z-components of the 
velocity and homogeneous suction for the y-component. In the spanwise and streamwise 
directions we employ periodic boundary conditions with periods L x and L z , respectively. 
The upper wall modifies the downstream velocity in the profile slightly. It now is 



u 



(lT(l -V S , 0), (2.5) 

where = E/<x>/(l - e h / s ). For h ->• oo the denominator is 1 and = Uoo. The 
99% boundary layer thickness slightly depends on h; it is <5 g9 o/ = 4.6055 for h — > oo and 
^99% — 4.6015 for h — 10, the box height used in most of this work. The laminar velocity 
profile uo(y) for this height is shown in figure [I] 

2.2. Numerics 

For all DNS (Direct Numerical Simulation) in this work we use the open source chan- 



nelflow library (Gibson 2012), developed and maintained by John F. Gibson. The code 
solves the full incompressible Navier-Stokes equations in channel geometries with peri- 
odic boundary conditions in the streamwise and spanwise directions and no-slip boundary 
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Figure 1. Laminar velocity profile uo(y) for a boundary layer with homogeneous suction Vs 

and for h — 108. 

conditions at the walls. It uses spectral discretization with Fourier x Chebyshev x Fourier 
modes in the x, y and z directions and finite differencing in time. The code has been 
modified to allow for the blowing/suction boundary conditions. As suggested by the chan- 
nelflow manual, computations in this work use 2/3-dealiasing in the x, z transforms and 
calculation of the nonlinear term in the rotational form (u-V)u = (V x u) x u+ V(u-u) /2. 
For time-stepping a third-order semi-implicit backwards-differentiation method is used. 

The resolution necessary to resolve all relevant flow structures has been carefully in- 
vestigated. In a box of size L x x h x L z = An x 10 x 2ir at Reynolds number 400 we 
increased the number of modes until all parameters of two edge states found at two 
different resolutions were identical to within 1%. A resolution of 48 x 129 x 48 modes 
was found to be sufficient. The high number of required Chebyshev modes is caused by 
the need to resolve small fluctuations near the lower wall. For wider and longer boxes 
the resolution was chosen such that the number of modes per unit length is at least the 
resolution demanded by the above analysis, i.e. M x /L x ^ 48/4-7T and M z /L z ^ 48/27T. 
We have examined the influence of the box height on the flow with a similar method and 
found that a height of h = 10 is sufficient to calculate edge states. It is not sufficient, 
though, for simulating fully turbulent flows, which need a much higher domain. 

We also found good agreement between our results and independent simulations by 
Schlatter and Duguet (2010, private communication). This way we could verify our im- 
plementation of both the ASBL modifications in channelflow and the edge state tracking 
algorithm. 



Edge state for a boundary layer 



5 



2.3. The edge of chaos and the edge state 

In the asymptotic suction boundary layer turbulent motion is observed while the laminar 
profile is still linearly stable. It is in that sense similar to other parallel shear flows like 
plane Couette flow ( Grossmann|20 00 ) or pipe flow (Eck hardt et aL|2007 ) and ideas about 
the state space structure can be carried over from these systems. There is a basin of at- 
traction around the stable laminar profile in which all initial conditions evolve towards 
the laminar flow. For sufficiently high Reynolds numbers and sufficiently large perturba- 
tions a second region appears in which the flow evolves towards a turbulent state (for 
our purposes here it does not matter whether this turbulent state is transient or persis- 
tent). In between those two regions is a laminar-turbulent boundary, the edge of chaos 
( Skufca et aL| |2006), which is characterized by initial conditions that neither decay nor 
become turbulent. The edge of chaos forms an invariant manifold of codimension one 
that separates the basin of attraction of the laminar state from the turbulent region of 
state-space. 

States within the edge of chaos evolve towards a relative attractor, the so-called edge 
state. This may be a fixed point as in plane Couette flow for low Reynolds numbers 
( Schneider et a/.|2008 ), a periodic orbit as in some parameter regions of plane Poisseuille 



flow (Toh & Itano 2003) and some low dimensional models (Skufca et al. 2006) or a 



relative chaotic attractor as in pipe flow ( Schneiderj|2007 Schneider et aL||2007 ), see also 



Vollmer et al. (2009). Its unstable manifold has codimension one, so that the edge of 



chaos is the stable manifold of the edge state. 

In order to find the edge of chaos and the edge state we use the technique described in 
Schneider eTal] ( |2007[ ) , see also |Toh fc Itano] ( |2003[ ) ; [Skufca et al.\ fl2006| ). For this edge 
tracking, we classify states that directly decay to the laminar state as lying on the "low 
side" of the edge and states that become turbulent as lying on the "high side". Starting 
with a turbulent state ut (on the high side) and the laminar fixed point uq, we use a 
bisection method to determine two states ul and ur on either side of the edge of chaos 
which are arbitrarily close to it. We can then propagate these states in time, always 
knowing that the edge lies between them. The two states quickly separate in the edge's 
unstable direction; after some time, we need a refining bisection to find a new pair of 
states that are closer together. The edge state tracking algorithm consists of a constant 
repetition of bisection and time integration, see figure [2] for a schematic representation. 

As indicated in figure [2j refining bisections are performed along the line connecting uh 
and uq rather than between uh and ul, because this method leads to better convergence. 
The reason for this improved performance seems to be connected with the folds in the 
edge, see |Moehlis et al\ ( |2004[ ); |Schneider eTai] ( |2007| ); | Vollmer et a!] ( |2009[ ): Bisecting 
between uh and the laminar fixed point seems to increase the probability of avoiding 
these folds and of staying near the edge. In our simulations we bring the relative L^- 
distance of ul and ur down to 10 -6 and perform a refining bisection once it is larger 
than f0~ 4 . 

We classify flow states as lying on the low or high side of the edge by thresholds of the 
cross-flow energy 

f 



E cf - 



dV, 



(2.6) 



L X L Z 8 , 

which measures the energy in the fluctuations perpendicular to the mainly downstream 
laminar flow. E c f should be a better indicator for the presence turbulence than just 
taking the L2-norm of the velocity field since it avoids the dominant contribution from 
the downstream component and focuses on the perpendicular parts that are necessary to 
sustain the three-dimensional turbulence. States with E c f greater than an appropriately 
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Figure 2. Schematic representation of the edge tracking algorithm. The black line represents 
the edge of chaos, the laminar fixed point uo is at the bottom. A bisection is performed on the 
line between the laminar and a known turbulent state. The states that bracket the edge are then 
propagated in time until their separation exceeds a threshold. Then, a refining bisection along 
the line connecting the new turbulent state and the laminar fixed point is carried out. 



chosen threshold become turbulent, while states with E c f smaller than another threshold 
directly decay. The use of E c f as a classification for states proved reliable and well-defined 
in our simulations, justifying the choice a posteriori. 



3. The edge state 

In this section, we discuss the results of edge state tracking at Re = 400 in a box 
of size L x x h x L z = 4n x 10 x 2tt and report the properties of the edge state. The 
resolution for this simulation was 48 x 129 x 48. A large number of edge state trackings 
for different initial conditions in this box have been performed. Some did not converge 
to any recognizable state (the dynamical properties of edge tracking will be part of a 
separate study), but all of those who did, converged to the same state, up to trivial 
translations - which we hence call the edge state. 

3.1. Description of the dynamics 

Figure 3]shows the energy in the in-plane velocity components, E c f, as a function of time. 
Only the trajectory eventually becoming turbulent is shown, the laminarizing one being 
visually indistinguishable. An outstanding feature of this plot are the periodic strong 
bursts during which the cross-flow energy E c f increases almost tenfold. These bursts 
provide the main starting point for understanding the edge state dynamics. The time 
interval between bursts is T = 1071 time units S/Uoq. 

Figure [4] shows the deviation of the flow from the laminar profile averaged over x at 
the minimum of E c f and exactly one burst later. To help visualizing structures across the 
periodic boundaries in z we show the box twice. This representation suggests that after 
T = 1071 the flow is the same except for a shift by n, corresponding to half a box width. 
After a second burst, the structures are hence back to their original location. Taking the 
changes in velocity into account, the full period is thus twice the time interval between 
bursts in E c t. 
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Figure 3. Time evolution of the cross-flow energy of the edge state. The initial transient be- 
haviour is not shown. The arrows at the top indicate the times when a refining bisection was 
performed, the circles indicate the times where the snapshots in the following figures have been 
taken. The locations of two minima are indicated with dashed lines; their distance is 1071 time 
units. 
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z z 

Figure 4. Fluctuating velocity field of the edge state at two successive minima of E c f. The 
velocity field is averaged in the downstream direction, the downstream component is colour 
coded, in-plane components are indicated with arrows with length proportional to the in-plane 
vector amplitude. The two images coincide when shifted by n = L z /2, indicating a periodic 
orbit with period 2142 time units. 



In order to study the dynamics of the edge state further, we present it in different 
representations at times corresponding to the circles in figure [3J Figure [5] shows snapshots 
of the fluctuating velocity field in a slice in the yz-plane. The plots are cut off at y = 2ir 
since in the upper part of the box no deviations from the laminar flow are visible. It 
becomes clear that the regions of high- and low-speed fluid already visible in the x- 
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Figure 5. The fluctuating velocity field in a slice in the yz-plane at x — at different times. 
Downstream velocity is colour coded, in-plane velocity is shown by arrows. The length of the 
arrows is proportional to the in-plane component of the velocity vector. The box has been cut 
off at y — 2%. The times correspond to the dots in figure [3] 



averaged plots are in fact streamwise high- and low-speed streaks. The snapshots hint 
at some internal dynamics, which below will be associated with a travelling-wave-like 
behaviour. 

The in-plane components in figure [4] suggest that vortices play an important role in 
the edge state's dynamics; we visualize them using the A2 vortex detection criterion 
(Jeong & Hussain 1995). Figure [6] shows the position of vortices by isocontours of A2 
and their interaction with the low-speed streak (isocontour u = —0.33, light gray) in a 
parallel projection, seen from the top of the box. Since the strength of the vortices varies 
considerably, we have normalized A 2 at every time slice such that min{\2) — —1 and 
plot the isocontour A2 = —0.6. It is coloured according to the sign of the downstream 
vorticity lu x in order to distinguish vortices that rotate left and right. 

We start describing the sequence of events at time t = where E c t is minimal. The 
flow consists of two streamwise streaks: one low-speed streak that reaches out to approxi- 
mately #99% into the flow, and beneath it one smaller high-speed streak closer to the wall. 
They are mostly aligned in the downstream direction, but the downstream modulation 
in figure [6] shows that it corresponds to a travelling- wave- like structure. The two streaks 
are accompanied by two streamwise elongated counter-rotating vortices, located at both 
sides of the low-speed streak and at different streamwise positions, just like the typi- 
cal flow structures identified in Jeong et al. (1997). The vortices create and sustain the 
streaks by linear advection of the laminar base flow. On one side a vortex lifts low-speed 
fluid from near the wall upwards and supports the low-speed streak. On the other side, 
fast moving fluid from the upper flow regions is pushed towards the wall where the base 
flow is slower, thereby creating a streak of high-speed fluid. Since the two vortices are 
counter-rotating, they both act to feed both streaks. For a long time the streaks persist 
without much change in shape while the energy in the in-plane velocity components E c f 
remains small. The vortices slowly grow in strength and start to "lean" over the low- 
speed streak, which in turn becomes wavier and the oscillations increase considerably. 
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Figure 6. Snapshots of vortex-streak-interactions in a view from the top. We show the isocon- 
tour u = —0.3 3 of the deviation from the laminar profile in light gray together with isocontours 
of A2 = —0.6 ( Jeong fc Hussain||1995[ ), where we have normalized A2 such that min(\2) = —1. 
The A2 isocontours are coloured red and blue for positive and negative downstream vorticity, 
respectively, with red indicating a vortex that rotates clockwise. Between t = 766 and t — 800 
the vortices switch places, thereafter new streaks are created. 
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Figure 7. Time traces of the edge state, (a) A plot of the edge state trajectory in a projection 
of state space. The circles correspond to the minima of E c f, crosses are drawn every 20 time 
units. The trajectory is a closed orbit passing through two minima. After a minimum of E c f 
the dynamics is rather slow, followed by a rapid motion as the next minimum is approached, 
(b) Time evolution of v at (0, 1.46, 0). The curve shows that two timescales are involved in the 
process. The inlet shows the same velocity component at time t and t + 2142 and reveals that 
the orbit is a periodic orbit without shift in the downstream direction. 

Approximately 700 time units after the burst the vortices lie on top of the now very con- 
torted low-speed streak. A little later the oscillations become too strong and the streaks 
break up - the event corresponding to the strong burst in E c f. Figure [6^d) suggests that 
the vortices actually switch places. The low-speed streak does not vanish in the process, 
but rather is torn apart into two smaller streaks, which then reunite across the periodic 
boundary conditions. After the break-up, the pairs of streaks and vortices are located 
at a spanwise location that is shifted by exactly L z /2 or half a box width. Comparing 
figures [6]ja) and (h) one finds that there is no shift in the downstream position of the 
structures. 

Further evidence for the edge state being a periodic orbit is presented in figure [7| 
The reduced state-space plot of two expansion coefficients in figure [7]Ja) shows that the 
trajectory is a closed orbit in this representation. The circle in the figure indicates the 
point of lowest E c f . The spacing between the crosses that are drawn at equidistant time 
intervals of 20 units shows that the trajectory moves rather slowly near the points, before 
rapidly going over to the other side in an event corresponding to the burst. In figure mb) 




we show the time evolution of a single in-plane velocity component at a fixed position 
in the box, namely v(0, 1.46, 0). The bursts show up as a slow drift from positive to 
negative values and back. An interesting feature are the short period oscillations. By 
averaging over 20 maxima of u(0, 1.46, 0) we find that the period is not constant during a 
cycle: directly after the burst, it is 18.4, corresponding to a phase-speed of approximately 
2/3t/oo. The period increases slightly and just before the burst it is 20.6. The inset in 
figure ^b) contains a magnification of the indicated black rectangle. In addition to the 
original curve, the dashed line shows the component at time t + 2142: that it lies exactly 
on top of the first curve demonstrates that the edge state is a periodic orbit and that 
there is no phase-shift - the high frequency oscillations are perfectly aligned so that 
the fast and slow dynamics are phase-locked. The velocity field during one of the short 
periods is reminiscent of the travelling waves found also in other systems. 

A space-time plot of the streak positions shows that the shift distance is always exactly 
L z /2. To determine the instantaneous position of the streaks, we average the fluctuating 
downstream velocity component u over x and y and take extrema as points of reference for 
the location of the streak. The velocities are fitted with a parabola to locate the minima 
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Figure 8. Localization and dynamics of low speed streaks, (a) The velocity profile of the edge 
state averaged over x and y at time t = (circles) and t — 850 (triangles). The position of 
the streaks is determined by fitting a parabola through the minima of those curves, (b) The 
z-position of the low speed streaks as a function of time. During a burst, the low speed streak 
splits into two streaks on either side, which then reconnect to form a new one. 



with higher precision. Two examples of the averaged profiles are shown in figure ^ a); 
the circles show the profile at t = 0, the triangles at t = 850 during the burst. While 
in the former curve, only one low speed streak is present, the latter one shows two low 
speed regions. In figure [8^b) we show the positions of the low-speed streaks as a function 
of time. That provide a reliable method to determine the shift distance of the structures 
at a burst and confirms that this distance is exactly it = L z /2. Moreover, the plot shows 
that the streak does not shift to one side at a burst but that two streaks form to the left 
and right of it which then reconnect to form a single one at the new position. 

The edge state obeys a shift-and-reflect symmetry, that is it is unchanged by a down- 
stream shift of L x /2 combined with spanwise reflection at an appropriately chosen 
plane. Formally, it is invariant under the operation cr[u, v, w](x, y, z) — [u, v, — w](x + 
L z /2, y, —z). We emphasize that no symmetry has been enforced in the calculations, but 
that edge trajectories are drawn to this apparently attracting invariant subspace. It is 
actually this symmetry that fixes the distance of the shifts: the only other possible lo- 
cation for the structures in the same invariant subspace is the one related by a shift of 
L z /2. 



3.2. A regenerating cycle? 



Hamilton et al. (1995) and Waleffe (1995) discuss a three step regeneration cycle for 



near wall turbulence. This quasi-cyclic process is based on the interaction between two 
alternating streaks and two streamwise elongated vortices. The streamwise vortices create 
the streaks by linear advection of the laminar base flow. But the streaks are linearly 
unstable to wavy modulations, that eventually lead to their break-up, causing a peak 
in the turbulent energy production. The resulting flow has a strong ^-dependence and 



nonlinear interactions lead to the recreation of streamwise vortices. Moehlis et al. (2004) 



were able to capture the essential features of the process with a nine-dimensional model. 
A similar mechanism has also been used to describe a self-sustaining process for coherent 
structures such as travelling waves flWaleffe|[l997l ). 

In the edge state in the asymptotic suction boundary layer, we see structures that 



closely resemble the travelling waves found in pipe flow ( Faisst & Eckhardt 2003| Wedin & 
Kerswell|2004)|Hof et a/.|2004[ ) or plane Couette flow ( |Nagata|1997||Clever fc Busse|1997 



Wang et al. 2007). They have rather high-frequency oscillations, see figure |7jb) . Such 
coherent structures have often been explained with the self-sustaining process (IWaleffe 
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Figure 9. Modal energies M(0, /3) and M(a, 0) as a function of time for a full cycle of the edge 
state. The changes in the two quantities are anticorrelated during most of the process, except 
for a small time interval t = 850 to 880. There are two phases, a long one during which M(0, /3) 
decreases and a short one where it increases. 



1997), but in the edge state there is a second time scale from the bursts. We find that the 



qualitative description of the low-frequency process of the bursts is more reminiscent of 



the process described in Hamilton et al. ( 1995 ) than the high-frequency travelling wave 



part. Important quantities for identifying the cycle are the modal energies 



M{k x 



ma. 



1 -J 1/2 

[u 2 (ma, y, n/3) + v 2 (ma, y, n(3) + w 2 (ma, y, n/3)] dy 



(3.1) 

The general interpretation is that M(0,/3) corresponds to the streaks and M(a,0) to 



the waviness of the streaks. Hamilton et al.\ ( 1995 ) found that during the cycle M{a, 0) 



and M(0, /3) are strictly anticorrelated. In their regeneration cycle two phases of ap- 
proximately the same duration are clearly distinguishable, one where the streaks grow, 
dM(0,/3)/di > 0, and one where streaks break down, dAf(0, (3)/dt < 0. We present the 
corresponding plot of the two quantities in figure [9] Af(0,/3) and M(a, 0) are anticorre- 
lated for most of the time, with the exception of a very small time interval during bursts 
- possibly indicating that the mechanism that leads to the recreation of the streamwise 
vortices is a different one. The streak formation phase after a burst takes only 200 time 
units, their slow breakdown lasts for more than 800 time units. In the discussion of fig- 
ure [6] in the previous section, we argued that streak formation is due to linear advection 
by the streamwise vortices. A sinuous streak instability that leads to a waviness of ini- 



tially straight streaks has often been observed, see for example Alfredsson et al. (2001) 



and there is little doubt that this linear instability is present in our case. Finally, the 
interactions leading to the recreation of the streamwise vortices are unquestionable non- 
linear - which means that all ingredients of the regeneration mechanism described by 



Hamilton et al. ( 1995 ) are present. 



3.3. The edge state in state space 

The dominance of the streaks, the long period and the dramatic bursts are reminiscent 
of behaviour seen in other boundary layers ( Kline et aL||1967 l 

In a low dimensional model obtained by proper orthogonal decomposition, |Aubry et al.\ 
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(1988) identified an attracting heteroclinic cycle which connects two fixed points, that 
are related by a discrete symmetry. The fixed points have exactly one unstable eigen- 
value whose real part is smaller than the modulus of the real part of the least stable 
eigenvalue. Trajectories spend a long time near the fixed points, before following the het- 
eroclinic connection to the other fixed point in a wild burst. Motivated by these studies, 



Armbruster et al. ( 1988 ) undertook further analysis of the phenomenon and concluded 



"that the heteroclinic cycles present in this model arise as a natural feature in the context 
of evolution equations which are translation and reflection invariant with respect to a 
spatial direction." ( Armbruster et al]|1988 p. 257) Since the required 0(2) symmetry is 
also present in our system due to the periodic boundary conditions in the spanwise direc- 
tion, we tried to connect our observations to their mechanisms. But despite considerable 
effort to locate the fixed points visited transiently by the trajectory, we have not been 
able to find them with a Newton search \ Viswanath|2007 ) . It eventually turned out that 
the mechanism responsible for the bursts is different, as we now explain. 

The origin of the slow oscillations was found by embedding the flow in a family of 
flows, obtained by keeping [Too fixed and varying Vs, thus exploiting the connection from 
ASBL to plane Couette flow (pCf ) . We work in a box of size Aw x 8 x 2tt at a resolution 
of 32 x 65 x 32 in order to reduce the computational effort. This resolution is a little 
lower than for the calculations above, but the edge state dynamics is still reproduced, 
though the values of some variables may vary slightly; for example the bursting period 
T increases by less than 1%. Since the edge states in both the ASBL as well as in pCf 
have a shift-and-reflect symmetry, we enforce this symmetry for all calculations in this 
section. 

We start in a plane Couette system with the lower plate at rest and the upper plate 
moving with velocity Uoo, where the Reynolds number based on half the velocity dif- 
ference and half the gap width is Re p c/ = UooH/Av = 800. The edge state is a fixed 
point in an appropriately chosen frame of reference and hence a travelling wave in the 
lab frame. The eigenvalues of the equilibrium state are calculated by means of Arnoldi 
iteration (implemented in channelflow by John F. Gibson); it has exactly one unstable 
direction, i.e. it is an attractor for edge trajectories. 

The suction velocity Vs serves as a homotopy parameter between pCf and the ASBL. 
As Vs increases from 0, the laminar profile gradually transforms from a straight line 
to an exponential, figure 10 a), while the boundary layer thickness 6 — vjVs decreases 
from oo to 1. Starting with the relative equilibrium that is the edge state in the plane 
Couette system, we continue the solution in Vs; it vanishes in what looks locally like a 
saddle-node bifurcation at Vs/Uoc = 6.418 • 10~ 4 . Restricted to the edge the edge state 
is a node and it corresponds to the upper branch of the saddle-node bifurcation. We 
were able to continue it to the lower branch which has two unstable eigenvalues. The 
bifurcation diagram E c f versus Vs is shown in figure 10 h). 

In figure [TT] we present the results of edge state tracking for four values of Vs after the 
critical value. Directly after the bifurcation, periodic bursts set in with an extremely long 
period. The period rapidly decreases with increasing distance from the critical suction 
value, as can be seen in fig ure [PT^b-d). A t every burst, the system shifts by L z /2. 



In view of the results of Aubry et al.\ ( 1988 ) and Armbruster et al. (1988) as well as 



Halcrow (20081 and Gibson et al. (2008&) it is a valid speculation that the edge state 



before the bifurcation is heteroclinically connected with its cousin that is shifted by L z /2 
and that the connection is attracting. Since both states are attractors (inside the edge) , a 
saddle must lie between them on this connection. For symmetry reasons, the heteroclinic 
connection exists as well in the other direction, making it an attracting invariant cycle. 
The resulting state-space structure is shown in figure |l2"](a); the unstable direction of the 
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Figure 10. Connecting the asymptotic suction boundary layer to plane Couette flow, (a) The 
laminar velocity profile uo(y) for different values of Vs as indicated in the legend in a box of 
height 8. For plane Couette flow (Vs — 0) the laminar profile is a straight line. It is smoothly 
transformed into the exponential profile of the ASBL as Vs increases, (b) Bifurcation diagram 
of the edge state (full circles). At Vs = 6.418 • 10 -4 , corresponding to S = 3.89, it collides with 
a lower branch equilibrium (open circles) in what looks locally like a saddle-node bifurcation. 
Restricted to the edge, the upper branch equilibrium is stable while the lower branch has one 
unstable eigenvalue. 



edge points perpendicular to the plane. The global bifurcation is hence a saddle-node 
infinite-period (SNIPER) bifurcation ( Strogatz||1994 Abshagen et a/.||2005 ), also known 
as saddle-node-on-invariant-cycle (SNIC), where the period of the cycle diverges near the 
bifurcation point (figure 12 b,c). The strong increase indicative of a divergence of the 
period can be read off from figure [TT] 

Further evidence for this scenario comes from a study of the explicit identification 
of the heteroclinic connection. The procedure we used to identify the connection is as 
follows: restricted to the edge of chaos, the node is an attractor while the saddle has one 



unstable direction. According to figure 12 an edge-trajectory that starts near the saddle 
will follow the heteroclinic connection either to the node or the shifted node. Using the 
edge state tracking algorithm to restrict the trajectories to the edge, we start from initial 
conditions near the saddle. A projection of state-space in figure [13] shows the two pairs of 
node and saddle, the heteroclinic connection to both sides as well as an edge state after 
the bifurcation. 

The long phase between two bursts of the edge state in the ASBL does hence not 
correspond to the transient visit of a weakly unstable fixed point but rather to the ghost 
left behind by a saddle-node bifurcation on an invariant cycle. 



4. Variation of flow parameters 

We now turn to a study of the properties of the edge state at different Reynolds 
numbers and in different domains. While there are some changes in detail, the edge 



Edge state for a boundary layer 
(a) V s = 7 ■ 1(T 4 (b) V s = 9 • 1(T 4 



15 




2 4 6 8 10 12 
T/10 3 
(c) V s = 12- 10- 4 



i 1 1 1 r 



J 



2 4 6 8 10 12 
T/10 3 
(d) V s = 15 • 10- 4 



i 

o 





i r 


1 — 


1 


— i 




J 


J 







4 6 8 10 12 
T/10 3 




i 1 1 1 r 



1 Ij Li —i 



4 6 8 10 12 
T/10 3 



Figure 11. Cross-flow energy of edge states for several values Vs- After the bifurcation point, 
periodic bursts set in with an extremely long period. The period decreases rapidly as the distance 
from the bifurcation point increases. 




c) 



Figure 12. Schematic representation of the bifurcation in state-space, a) for small Vs two pairs 
of fixed points, one stable (full circles), one unstable (open circles), exist which are related by a 
discrete symmetry, namely shifting by L z /2; they are heteroclinically connected as indicated, b) 
as Vs increases, the fixed points annihilate in a saddle-node infinite-period (SNIPER) bifurcation 
c) The fixed points are gone for larger Vs, the edge state now goes around the cycle. For Vs 
near the bifurcation value, the period of the cycle diverges due to the ghost left behind by the 
saddle-node bifurcation. 



states' structure is preserved in all cases. It consists of two streamwise streaks and two 
streamwise elongated counter-rotating vortices and there are periodic bursts where the 
state is shifted by half a wavelength. 

The criterion for the convergence of a calculation is based on the cross-flow energy. We 
only consider states where the evolution of E c f is in some sense recognizable, for example 
constant or periodic. We found that all states falling into this category display periodic 
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0.06 



Figure 13. The heteroclinic connection in state-space at Vs = 5 • 10~ 4 in a two dimensional 
projection. We show the node (full circle) and the saddle (open circle) together with the pair of 
fixed points obtained by a shift of L z /2 (the projection is independent of x). The black lines are 
obtained by edge tracking from initial conditions which are slight perturbations of the saddle 
point; they reveal the heteroclinic connections from the saddle to both stable fixed points. The 
dotted line is the (bursting) edge state at Vs = 7 • 10~ 4 , slightly above the bifurcation. 



bursts. Convergence is assumed if E c f has at least three peaks where their height and 
the time between the peaks agree within a relative tolerance of 2%. 

4.1. Reynolds number 

In this section the Reynolds numbers is varied while the box size is kept fixed at 47r x 2tt. 
Calculations still take place in the transitional regime, with Reynolds numbers between 
320 and 750. Below Re = 320 the cross-flow energy at bursts is of the same order as 
the cross-flow energy of fully turbulent states and the distinction between the edge state 
and turbulent motion becomes difficult. The highest Reynolds number we calculated is 
750. There is, however, no indication that the structure of edge states would change for 
slightly higher Reynolds numbers. 
Figure 



14 



shows the cross-flow energy E c t as a function of Reynolds number. The 
maximum value of E c f at a burst is plotted with triangles. This value may be taken as 
a measure for the burst intensity, it decreases monotonically with Reynolds number. A 
potential explanation comes from the observation, that at higher Reynolds numbers the 
critical amplitude for triggering transition to turbulence is lower ( Levin et aZ.||2005 ) and 
that the edge state's energy reflects this descending threshold. The circles show the time 
averaged cross-flow energy. It decreases monotonically and smoothly with Re for the same 
reasons as the maximum energy. A power law gives the best fit with (E c f) T ~ i?e~ 3 , 
plotted as a dashed line in the same figure. 

Figure 15 'a) shows the minimum of E c f which corresponds to the undershoot after a 
burst. As in the former figure, E c f . mostly decreases with Re. In 15 b), we present 
the bursting period, the time between two bursts. It increases almost linearly with Re for 
most Reynolds numbers. The states spend longer times between two bursts, the bursting 
process itself does not take longer. 
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Figure 14. The energy of the cross-flow E c f as a function of Re. Green triangles: maximum of 
E c f at bursts; red circles: E c f averaged over one period. The dashed line is a fit to a power law, 
E c f ~ Re . The energy generally decreases with increasing Re. 
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Figure 15. Variations with Reynolds number, (a) The minimum cross-flow energy as a function 
of Reynolds number, (b) The bursting period. Both plots show a discontinuity at Re = 425 and 
at Re — 455, marked with dotted lines. Apart from that, the variations are rather smooth. 



We conclude this section by reporting a feature in the range of Reynolds numbers 
between 425 and 570 that awaits explanation. At Re = 425 the minimal cross-flow energy 
E c f min undergoes a drastic and sudden drop; after that it recovers continuously until 
Re = 570, where it again turns downwards. The bursting period rises discontinuously at 
Re = 425 and drops again at 570. We emphasize that the general edge state structure 
does not change: the bursts are periodic, the states shift by half a box width and there 
are no noticeable changes seen in the flow visualizations. 

4.2. Streamwise and spanwise domain-size 

In this section, we study edge states at Reynolds number 400 but in boxes of different 
size. L x varies from 3.27T to 77r in steps of 0.27T, with L z ranging from L x /2 — 0.57T to 
L x /2 + 0.4-7T in steps of 0.17T. In smaller boxes it was difficult to find sustained turbulent 
states. 

We present the results as pseudocolour plots of the quantity of interest, with L x on the 
horizontal axis and L z on the vertical one. The result for the box L x x L z a represented in 
the rectangle (L x ±0.1) x (L z ± 0.05). Rectangles where no calculations were performed 
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Figure 16. The time-averaged cross-flow energy for different box sizes. Overall the cross-flow 
energy is smaller in larger boxes. Boxes where no calculations have been performed or where we 
found no periodic bursts are left white. The dashed line indicates L x = 2L Z . 



or showed no periodic bursts are left white. The "diagonal" L x = 2L Z is indicated by the 
dashed black line. 



Figure 16 shows the time-averaged cross-flow energy. As a general trend, the cross 
flow energy is high in the lower left corner of the plot where the box size is small and 
low in the upper right corner, in larger boxes. A notable exception are 13 states in a 
small region around 5.47T x 2.57T. The large scale trend can be readily interpreted by the 
linkage of the flow strcutures to the boundary layer thickness. The streak size is limited 
by the boundary layer thickness, since no shear gradient required to sustain the streaks 
is present above. Therefore, the streaks' spanwise extent does not grow to fill the whole 
box. 

E c f is normalized by the boundary layer volume. A plot of the minimum cross-flow 



energy multiplied with the box area L X L Z in figure 17 shows that in contrast to figure 16 



no trend is present. To some extent, the minimal cross-flow energy is lower along the 
dashed black line L x = 2L Z . We think that this behaviour may be taken as a hint that 
naturally a wavelength ratio of 2 : 1 is preferred, as also seen in plane Couette flow 
dClever fc Busse||1997[ |WaJeffe|2003 ). 

The bursting period T, presented in figure 18 varies inversely to the cross-flow energy. 
It is short for small boxes and long for bigger boxes. There are probably at least two 
important contributions to this effect. The first one is given by the mechanism already 
explained in section |4.1| higher energies trigger the bursts more frequently. The anti- 
correlation is easily seen by comparing figures [16] and [18] But a second contribution is 
revealed by a closer inspection of figure 18 the bursting period along a horizontal line of 
constant box width is roughly constant, but it increases drastically along vertical lines, 
as the box width grows. In wider boxes, the oscillations of the streaks have more time to 
build up before breaking up by touching their neighbour. 

In figure 19 the cross-flow energy is plotted as a function of the boundary layer volume 
L x ■ L z . We distinguish three qualitatively different regions: small boxes with L X L Z < 
107T 2 , larger boxes with L X L Z > 10ir 2 and the states around 5.47tx2.57t, which are marked 
with open symbols in figure [i~9] In larger boxes, the variation of the cross-flow energy is 
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Figure 17. Minimum cross-flow energy for different box sizes. The energies are not normalized 
by volume. States near the dashed black line L x = 2L Z seem to have a lower value of E c f . . 
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Figure 18. The bursting interval T for different box sizes. T increases to the upper right corner, 
the trend is the reverse of the cross-flow energy, cf. figure 16. Furthermore, T is roughly constant 
along horizontal lines and increases for wider boxes. The dashed line indicates L x = 2L Z . 



rather smooth. E c f max is roughly constant in those boxes: the burst intensity does not 
depend on the box size if the box is large enough. Both E c j min and (E c f) T decrease for 
increasing box size, but the variation is very small for the largest boxes. In the small 
boxes, the minimum and time-averaged cross-flow energy are not constant as a function 
of box size. States in these boxes seem to depend heavily on the available space and 
the greatest possible wavelength. Many states in these boxes have a comparable value of 
E c f max which is roughly 50% higher than in larger boxes. 

The states from boxes around 5.4-7T x 2.57T, marked with open symbols, lie outside 
the rather smooth variation of the curves, both E c f max and (E c f)T are considerably 
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Figure 19. The cross flow energy E c f as a function of the boundary layer volume. Three 
qualitatively different regions can be distinguished: states with L x ■ L z < 10 are somewhat 
higher in energy and show wide variations. For larger values the maximum of E c f seems to 
prefer a value near 25 and the variations in all energies are reduced. The third group of states, 
marked with open symbols, are from the boxes around 5A-k x 2.5ir. They lie outside the general 
trend and seem to form a group of their own. Theses states include, among others, the period 
doubled state shown in figure [22] 



higher and also E„t . varies more. The burst intensity, E r t , in these boxes has 

° w mzn J ' W max 1 

approximately the same value as in the small boxes. There presumably exist two periodic 
orbits in these boxes, one of which is stable in small boxes while the second one is 
stable in larger boxes, where in turn the first one is weakly unstable. We assume for 
the open symbols, that the edge state tracking algorithm approaches the unstable first 
one and integration time is too short to see the weak instability and converge to the 
"true" edge state. As an example, in figure [20] we present an edge state tracking in 
the box 4.27T x 2ir where that happens. We first see 5 peaks of the same height and at 
equal distance, where the edge state tracking seems to be converged. But then the curve 
changes and another state is approached, where E c f and E c f are lower and higher, 
respectively. Figure [20] shows that sometimes edge tracking has to be continued for very 
long times before converged results are obtained. It would be possible to distinguish 
transient states from true periodic states if they could be identified numerically exactly, 
say using Newtons method. Then the differences in the number of unstable eigenvalues 
would show up. However, it has not been possible to achieve this because of the long 
periods. Given the large number of calculated states we feel that the overall trends 
discussed in this section are reliable. 
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Figure 20. Non-converged edgetracking evidenced by the cross-flow energy in a box of size 
4.2-7T x 2n. Up to a time of about 6000 the tracking seems to converge rapidly. However, beyond 
that the curve changes: the bursts are less intense, the bursting period is a little lower - the 
edge state tracking has converged to another state, even though the first one would have been 
considered as converged by the criterion described at the beginning of this section. 



4.3. Along the line L x = 2L Z 



A box size close to L x = 2L Z has been found to be optimal in plane Couette flow ( Clever 
&; Busse||l997 Wang et aL|2007 ) and also seems to be a good domain here. We therefore 
present a more detailed study along the diagonal L x 



2L, 



In figure 21 the cross-flow energy is shown. For calculations marked with full symbols, 
we have imposed the shift-and-reflect symmetry, for calculations marked with open sym- 
bols, no symmetry has been enforced. The overall trend is the same already discussed 
for the colour plots in the previous section: with growing box size the cross-flow energy 
content per unit volume decreases. For L x > 5tt, the minimum cross-flow energy settles 
at a low value. 

The higher resolution possible for this one dimensional cut through the space of length 
reveals non-mononities near L x /ir — 3.8, 4.44, 5.0, 5.72 and 6.4. Both E c f max and (E c f)T 
show a higher value than would be obtained by an interpolation, while E c f is slightly 
lower. The rather regular spacing between those points of about 0.647T to 0.727T hints at 
some kind of resonance phenomenon, but we could not yet pin this down. 

4.4. A period doubling 



Figure [22] shows the cross-flow energy of the edge state in the box 5.67T x 37r. The burst 
intensity, measured by the maximum of E c f, is no longer the same for every burst, it 
alternates between a lower and a higher value. The periodicity of E c f is not completely 
lost, but the time interval after which it repeats has doubled. The period of the periodic 
orbit that is the edge state is not affected since it is twice the bursting period. The edge 
states' structure is still preserved, at a burst the streaks are shifted by exactly L z /2. 



5. Conclusions 

The periodic orbit in the asymptotic suction boundary layer which we have character- 
ized here shows dynamics on two different time scales: similar to invariant states in plane 
Couette and pipe flow, it shows the typical advection of downstream vortices on a fast 
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Figure 21. Cross-flow energy at 5ft = 400 for domains with L x = 2L Z . Upper panel: maximum of 
E c f at bursts; middle: E c f averaged over one period; bottom: minimal value of E c f. Calculations 
where a shift-and-reflect symmetry was imposed are marked with full symbols. 
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Figure 22. Cross-flow energy of the edge state in the box 5.6n x Stt. The burst intensity is no 
longer the same for every burst, but alternates between a lower and a higher value. The dashed 
lines correspond to the two maxima and confirm the periodicity of the state. 



time scale. Regular energetic burst events, not unlike those that have been reported as a 
key feature of boundary layer flows, occur on a much slower time scale. Although the flow 
can be turbulent and the dynamics is characterized by positive Lyapunov exponents, the 
two time scales remain phase-locked over the whole period of several thousand advective 
time units. In view of the translational symmetry of the asymptotic suction boundary 
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layer these phenomena could well be the periodically continued version of the evolving 



boundary layer edge state described by Duguet et al. (2012[) 



For the bursting dynamics, an explanation from two different view points is given. 
Physically, the bursting is the result of vortex-streak interactions. A pair of downstream 
vortices creates and sustains the streaks by linear advection. The low-speed streaks are 
linearly unstable and develop a waviness, which in turn causes the vortices to become 
tilted with respect to the x-axis. The two counter-rotating vortices "lean" over the low- 
speed streak, tearing it apart, and eventually switch sides. Afterwards, the streaks are 
recreated by the vortices and the process begins anew, only shifted by half a box width. 
From a dynamical systems point of view, we are able to associate the bursting with a 
SNIPER-bifurcation: in plane Couette flow, two symmetry-related pairs of saddle and 
node are heteroclinically connected. As the homotopy parameter, the suction velocity, 
is increased, each of the pairs collides in a local saddle-node bifurcation, leaving behind 
a periodic orbit along the former heteroclinic connections. The discussed periodic orbit 
provides one mechanism for energetic bursts, adding one more possibility to the general 



bursting phenomena discussed by Robinson (1991). It remains an open questions under 



which conditions these phenomena can be traced back to the instability of low-speed 
streamwise streaks as in the edge state solution. 

The obvious spatial localization in the edge calculations in the Blasius boundary layer 
( Cherubini et aL|2011 Duguet et aL|2012 ), as well as similar structures in confined flows, 



suggests the existence of spatially localized counterparts of the periodic orbit discussed 
here. We expect that further comparisons between these and other features found in 
the asymptotic suction boundary layer and the full Blasius boundary layer will help 
to identify and characterize shared generic features and to separate them from suction 
related influences. 
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